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Abstract 


The systematic investigations of the role of the virtual charm axial currents 
in the decay of i3-mesons to p' K and the spin structure of the nucleon are 
performed. We reduce the divergence of the virtual charm axial current to 
a specific gluon operator. Since this operator receives the main contribution 
from topologically nontrivial components of the QCD vacuum, we rederive 
the Diakonov-Petrov Effective Action, based on the instanton QCD vacuum 
model. This action is applied to the calculations of the coupling of r]' to 
the charm axial current and found — me) = —(12.3 ~ 18.4MeV), 

providing a possibility for the explanations of the recent experimental data 
on B rj' A-decay. Analogous calculations of the virtual charm content of 
the nucleon spin leads to Ac(p ~ me) = — (0.015 ~ 0.024), which is one 
order of magnitude smaller than the analogous contribution of the strange 
quark. 
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I. INTRODUCTION 


The processes sensitive to the OZI rule violations and the contributions of the nonva¬ 
lence components of hadrons become of great interest. In the present paper, we would like 
to discuss the role of the charm content of the light hadrons in the 5-meson decay to rj' 
and i^-mesons and in deep-inelastic scattering with polarized charged leptons on polarized 
nucleons. In both cases there are no charmed hadrons both in the initial and the hnal states, 
which means that a charm may give a contribution only through virtual processes. 

Decay of 5-mesons to rj' K 

Recently there is a great theoretical interest [3-13] on the new experimental data on the 
branching ratios of the decays of 5 — rj' K[l, 2]: 

5r(5^ ^ rj'K^) = (6.5l};^ ± 0.9) x 10-^ (1.1) 

5r(5° ^ r]'K^) = (4.7l|^ ± 0.9) x 10“® . (1.2) 

In the Standard Model, the Cabibbo favored b —> ccs elementary process may be followed 
by conversion of cc pair into r]' through gluons. The amplitude of this process is described 
by 

M=^K,Oi(V(p)|c7;.75c|0)(ir(g)|s7;.6|5(p + g)). (1.3) 

Here Gp is the weak coupling constant. Vet, V*g are Kobayashi-Masukawa matrix elements 
and ai = 0.25 is the phenomenological number obtained from a £t to experiments (see [3] 
for the references). The matrix element 

(h'(p)|c7/.75c|0) = (1.4) 

is non-zero due to the virtual cc —> gluons transitions. Certainly, this matrix element 
is suppressed by l/m^ factor. However, due to the presence of strong nonperturbative 
gluon helds in the QCD vacuum together with the Cabibbo favored b ^ c transition the 
suggested b ccs mechanism (1.3) may be expected to compete appreciably with the other 
mechanisms of the 5 — Kg' process [9,12]. 

If we assume the dominance of the virtual charm mechanism (1.3), the branching ratio 
is written in terms of as [3] 

f(c) 

If]' 

1 GeV 

Using the experimental data, (1.1), it is found ~ 140 MeV {"'exp”). 

This value perfectly coincides with the estimate of Halperin and Zhitnitsky [3]: 

= (50 - 180) MeV. (1.6) 

On the other hand, a recent phenomenological study placed another bound on f^f\ namely 

(c) fc) 9 

—65 MeV < < 15 MeV, with being consistent with zero by analyzing the Q'' 


Br{B Kg') ~ 3.92 ■ 10 


-3 
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evolution of the form factor [7], and more recently it was estimated from observed 
ratio of J/'^ decay to rj' and rjc the value of = —(6.3 ± 0.6)MeV [8]. Other similar 
estimation which leads to < 12MeV was made in [13]. Ah et ah considered the 

complete amplitude for the exclusive 5—meson decays, including r]'K channels, where they 
combined the contribution from the process b s{cc) —> s{gluons) —> sr]^ ^ with all the 
others [5,6]. Their estimations gave \ ^ 5.8 MeV[5] and = —3.1 (—2.3) MeV (for 

rric in the range 1.3 - 1.5 GeV)[6] in agreement with the analysis[7]. They stressed the 
importance of the sign of and found a theoretical branching ratio in the range 

Br{B r]'K) = (2 - 4) x 10■^ 

which is somewhat smaller than the experimental one (1.1). The similar analysis made in 
[10] led the authors to conclude that = —50MeV might provide the explanation of the 
data. 

Having this situation, it is important to recalculate to clarify the mechanism of 
B —>■ Kf]' decay in the similar framework performed by Halperin and Zhitnitsky [3]. 

Polarized PIS 

In the constituent quark model, the spin of the proton is supposed to be carried by g = 
u, d valence quarks so that AS = Au + Ad = 1. The quark spin Aq is dehned as 

Aq2mNS^ = {p,s\qj^'y5q\p,s), (1.7) 

where ttin and are the mass and the spin of the nucleon, respectively. On the other hand, 

Aq = [ dx {qnix) - qL^x )), 

Jo 

where qR(L){x) are quark distributions of chirally right-handed (left-handed) quarks in a 
polarized proton. 

Deep inelastic scattering (DIS) with polarized charged leptons on polarized targets pro¬ 
vides the investigations of the quark distributions g_R(L). These quantities are extracted from 
the structure function gi{x, measured in polarized DIS using the parton model relation 
gi{x,Q^) = {qR{x) - qiix)). 

So called ’’spin crisis” problem is related with a very large disagreement between the 
experiments and the prediction of the naive constituent quark model for the hrst moment 
of the proton(neutron) spin structure function T^*'”^ dehned by 



Here, 

g\ = Au — Ad = 1.25 

and 

g\ = {Au + Ad — 2As) = 0.69 

are the isovector and the octet axial charges measured from neutron and hyperon decays 
and 

= AS = (An -|- Ad + As). 
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The expectation was g\ = 1 but the measured value of this quantity given by EMC in 1988 
implied that g\ ~ 0.12. It means that only ~ 12% of the spin of the proton is carried by its 
quarks! The modern value is g\ ~ 0.3 [15-17]. 

In all of the consideration of the nucleon spin the contribution of the charm was neglected 
completely. With the account of this contribution 

= AS + 2Ac 

and we have now the problem to estimate also Ac. We may expect sizable value of Ac, 
since in this case we are dealing again with axial currents which may be strongly affected 
by the vacuum nonperturbative gluon helds. The previous calculations of Ac by Halperin 
and Zhitnitsky [18] and Blotz and Shuryak [19] gave quit different results Ac ~ 0.3 [18] and 
Ac/AS = — (0.2 ~ 0.08) [19]. They contradicted each other in the sign and also absolute 
value. This is our motivation for the recalculation of the charm contribution to the spin of 
the nucleon. 

Axial currents of virtual charm 

The symmetry of the classical lagrangian may be destroyed by quantum fluctuations 
[20-22]. In gauge theories the axial anomaly arises from noninvariance of the fermionic 
measure against axial transformations in the path integrals of the theory[23] (see also ref. [24], 
concerning higher-loop corrections). The present problem is intimately related with this 
phenomenon. 

In the following we will work only with the Euclidean QCD, where its convention is 
written in the footnote^. In the Euclidean QCD the axial anomaly in the light quark axial 
current in the chiral limit reads 

where tfjf is the light quark field (/ = m, d, s) and g is the QCD coupling constant. 2GG = 
where Gl^ is the gluon held strength operator with a being the color index. 

The situation with heavy quarks is very different, since we must take into account the 
contribution of the mass term. The divergence of the axial current of charmed quarks has a 
form: 


2mcC^5C, ( 1 . 10 ) 

The hrst term in (1.10) again comes from noninvariance of the fermionic measure (or in 
other words - from Pauli-Villars regularization). The main problem here is to calculate the 
contribution from the second term in (1.10). It is clear that this one is reduced to the 
problem of the calculation of the vacuum expectation value of the operator 2 mccl 75 C in the 
presence of a gluon helds. 


^ixMo — xea, xmi — XEi, Amo — iAE4, Ami — —Aei, V'm — V’e, — ^e, 7mo — Iea 
iMi = ilEi: iMb = lEb- the following we will omit index E. 
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In the path integral approach, the calculation of the contribution of this term to any 
matrix element over light hadrons may be considered in the sequence of the integrations. In 
the first step, the integration over c-quark is performed, and the next step is the calculation 
of the integral over gluon helds and dually the integration over light quarks. 

We consider here the first step - the integration over c-quarks. The resulting (2mcC^75c) 
must be a gauge invariant function of the gauge held A and therefore must be expressed 
through the gluon held strength tensor and their covariant derivatives. The matrix elements 
of such types of the operators are almost completely defined by nonperturbative topologically 
nontrivial (like instantons) contributions at least for the small virtualities flowing through 
these operators, which is our case. This one can be clearly justified by the consideration of 
low-energy theorems for the various matrix elements of the gluonic operator in (1.9) [25]. 
It was shown that only with account of the instanton-like configurations of vacuum gluonic 
fields, it is possible to satisfy these theorems and that Diakonov & Petrov (DP) chiral quark 
Effective Action [28,29] based on the QCD instanton vacuum model [32,33] reproduce well 
these low-energy theorems in the chiral limit but fail beyond this limit [26,27]. To our 
present knowledge the instanton structure of QCD vacuum is characterized by the average 
size p and by the average inter-instanton distance R, which are [33, 32] 

p = 1/3 fm, R = 1 fm. (1.11) 

Therefore the packing parameter (p/i?)^ = 0.012 is small, legitimatizing independent aver¬ 
aging over positions and orientations of the instantons. 

In the present paper we first calculate gluonic operators in the divergency of the vir¬ 
tual charm axial currents, further we rederive the DP Effective Action starting from the 
Lee&Bardeen result for the quark propagator in the instanton media and apply this action 
for the calculations of the correlators with gluonic operators. Einally we apply these results 
to the calculations of the virtual charm effects in above mentioned problems of the 5-decay 
and DIS. 


II. THE DIVERGENCE OF THE VIRTUAL CHARM AXIAL CURRENTS 

We calculate the expectation value of the operator 2mcC^75C in the presence of gluon 
field. Eor this, we take the method, which was developed by Schwinger in electrodynamics 
many years ago [20], and later was applied to QCD by Vainshtein et al. [34]. The key point 
of this method is based on an assumption of a possibility of an expansion of the Green 
function such as ( 2 mcC'l’ 75 c) over G/m^. 

We introduce first the coordinate and momentum operators, and respectively, 
which satisfy [p^,Xy\ = , [p^,Pu\ = [X^,X^] = 0. We dehne then the covariant 

momentum operator satisfying the following commutation relations, 

[V,, X,] = , [V,, VA = igGlX , ( 2 . 1 ) 

where R is a generator of the color group and is the gluon field strength tensor. More¬ 
over, we introduce a formal complete set of states |a;) as the eigenstates of the coordinate 
operator 

XX) = ^X) > ( 2 . 2 ) 
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which satisfies 


{y\x) = — y) , j d'^x\x){x\ = 1 . (2.3) 

In this basis, the operator V^J, acts as a covariant derivative 

{y\V^,\x) = iDf,{y\x) = (id^ + gAl{x)t^) 5^^\x - y) . (2.4) 

The algebra (2.1) is the basic tool of the Schwinger formalism. We expand the Green 
fnnctions in the gluon background field, and need to use only this algebra in each order of 
expansion. 

Next, we calculate the expectation value of the operator 2mcC^{x)'-^^c{x), which is the 
Green function that should be expanded by using the Schwinger method. We consider the 
integration over c-quarks. In the path integral approach, we define : 

{2mcc\x)'y^c{x)) = J ZlcZlc^ 2 mcC^(a;) 75 c(a;) exp |y d'^y c\y){V + imc)c{y)'^ (2.5) 

with V = Since the argument of the exponential is quadratic in the quark field c, the 

path integral (2.5) can be written in the form, 

(2mcc'^(a;)75c(a;)) = 2mcdet IIP + imdl (a;|Tr75 ^^—^ -|a;) , (2.6) 

V + irric 

where Tr denotes the trace over spin and color indices. Since a; is a continuous variable, the 
operator has an infinite number of matrix elements, and in calculating the determinant of 
this matrix there arise infinities of various types. Hence, the determinant, det ||P + irricW, 
must be regularized in the standard manner as 

(P + imc){p + iM) 
ip + imc){V + iM) 

where M is the Pauli-Villars regulator mass. Eq.(2.6) must be gauge invariant and expressed 
through the gluon field strength tensor and the covariant derivatives. 

Eq.(2.6) is expanded in the series of a power of G/m^ under the assumption that the 
field strength is much less than the square of c-quark mass Here, we will take into 
account 0(G‘^) and 0{G^) terms in the expansion of (2.6). We start from the calculation of 

H{x) = 2mc{x\TT'y5——^ -|a;) . (2.7) 

P + irric 

The calculations of the det ||P+imc|| is presented in the Appendix, since it give a contribution 
to (2.6) starting 0{G*) terms. 

Using the formulas 

p 2 _ -p2 _|_ j _ -- /p _ \ 

2 V-irric 

where aG = = |[7/i,7i/] and / is an identity matrix, Eq.(2.7) reduces to 
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1 


1 


H{x) = 2 mc(a;|Tr 75 - 


V + inic V — irric 
= - 2 im 2 (a;|Tr 75 -—- ^ „ J x) 


{V — imc)\x) 


'V^ + ml + ^aG' 


= —2iml{x\Ti'y^ 


9 


aG- 


:aG- 


+ ml 2 + ml 2 + ml 


--aG- 


1 


:-(tG- 


1 


--aG- 


V‘^ + ml2 V^ + ml2 + ml2 + ml 
= H2{x) + Hs{x) + --- . 


+ 


lx) 


( 2 . 8 ) 


It has been used here that the trace of an odd product of 7 matrices vanishes. 

It is straightforward to calculate the second term of the right hand side of (2.8), if neglect 
the noncommutativity of the operators and G^^ in (2.8). Thus, 

= ^(-|Tr75(^y^(aG)3|x) . 

In that case, since the operator can be replaced by the ordinary momentum we can 
use the evident formulas. 


(X 


(■p^ -I- fYi'iy 


|x) 


1 

(27r)^ (P 2 + m^)” 
d‘^p 1 
(27r)^ (p^ + ml)"- 


A=0 


2'^'K‘^{n 


1 

l)(n — 2)mc*'"'~^^ 


(2.9) 


and 


Tr75(aG)3 = 2^* tr.GGG = -2^fabcG^G^G^ , 


where GGG = G^uG^aGa^j. and fate is the structure constant of the color group. As a result. 

We still have to calculate H 2 {x) which contains the gluon held strength G^^ to the 
second power. However, the calculation of H 2 {x) needs much more efforts owing to the 
noncommutativity of the operators, which is not negligible. For the sake of the systematic 
momentum integration of H 2 {x), one must transfer all operators containing to the left 
hand side in the trace. Using the following expansion : 




+ ml 


+ ml 


P 2 + ml 
1 

+ 


aG 


aG 


p 2 + m 


[P^aG] 


p 2 + m 


:[P^aG] + 


(P 2 + 


2 )^ 


(p2+^2)2C > J (p2 + ^2)3 

[P^[P^[P^aG]]] + ■•• , 


[P^[P^aG]] 


( 2 . 11 ) 
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we can rewrite H 2 {x) as 


ig^ml 

H2{x) = -— (a;|Tr75 


i(jGY + 


{V^ + mlY 

= hi{x) + h 2 {x) + h'i{x) + 


(3lp2, l■p^ aG]]aG + 31V\ aGf + ffG|P^ |P^ ffG|]) + ■ ■ 'j |i> 

( 2 . 12 ) 


The calculation of the commutators is performed systematically. For an arbitrary operator 
Q, which satishes Q] = iD^Q, 

[V\Q] = D^Q + 2tV^D^Q . 

Repetitive use of this identity leads to 

[V^,aG] = D‘^aG + 2iV^,D^aG (2.13) 

[V^, [V^, aG]]= D^aG + 2iD^G^^ ■ D^aG + 2iV^{D^D^ + D^D^)aG 

—A'PyGi,^D^aG — AT’y'P^DyD^aG . (2-14) 


In the following, we assume the quasi-classical vacuum and the source p, which is dehned 
by \P^iG^y\ = iD^G^y = p, is sufficiently small. So, the second term in (2.14) may be 
negligible. Since other higher commutators contribute in H 2 {x) as at least third power of 
Gfj,ui we can neglect the corresponding terms which are denoted by dots in (2.12). 

The momentum integration of each term can be performed easily since we need only both 
0{G‘^) and 0{G^) terms. The calculation of the first term in H 2 {x), hi{x), is somewhat 
technical. By using the translational invariance of 

(a^|Tr75 (aG)^) 

with an arbitrary momentum and extracting all g^-terms as in Appendix, hi{x) is given 
as 

hAx) = + 0(G‘) , (2.15) 

However, 0{G‘^) term may be neglected for our purpose. So we may regard the calculation 
of hi{x) as an integral over an ordinary momentum p^, as we have calculated in (2.9). The 
term hi {x) contributes to cancel the axial anomaly from the noninvariance of the fermionic 
measure. 

The second term in H 2 {x), h 2 {x), is rewritten as 


h2{x) = 


• 2 2 
%g 


• 2 2 
%g 


T|Tr 75 


{x\ 


{V^ + m 
1 

{V^ + rrilY 


lY 


{2[V^, aG]aG + aG[V^, aG]) \x) 


\x)"Ii'y^aG ■ D^aG 


(2.16) 

(2.17) 


+6i {x\ 


(P 2 + mlY 


P^|a;)Tr 75 (TG ■ D^aG \ + (Total derivative) 



We can omit here also the terms which contain a single operator V^. The reason is that the 
matrix elements, {x\{V'^ + must be denoted in terms of G^i, and D^. The first 

nonvanishing term which can give a contribution to this matrix element is D^G^. It is clear 
that it leads to 0{G^), which we do not calculate here. 

By using the Bianchi identity, it is easy to show that 

D^G^u = —ig [Ga^i, Gau] + D^D^Gau — DyDaGan ■ (2-18) 

where the second and third terms in the right-hand side may be neglected. The reason is 

'~^^'^5^G(Jfii/Dfj^DciGau \p^fJLV ' ^^cGXpDpiDaGav 

av 

= 4:exppuGcDpGxp ■ D^Gau + (Total derivative) 

= StTcDpGpi, ■ DaGau + (Total derivative) 

= 0 -|- (Total derivative), 

because of the Bianchi identity, D^G^j, = 0. Then, D'^Gp^, is the second power of G, and 

D^aG = —‘2ig(jpyGapGciv ■ (2.19) 

The solution of h 2 {x) is obtained as 

h 2 (x) = = -H.Xx) . (2.20) 

Hence, h 2 {x) cancels with H^{x) in (2.8). Thus, it will be only the term h^{x) in H 2 {x), 
which contribute to the divergence of the axial current of charmed quarks. 

Our remaining work is to calculate only the third term in H 2 {x), namely, 

hsix) = (3[P^ [P^aG]]aG 

+3[iP^ aG]^ + (jG[V^, [V\ (tG]]) \x) . (2.21) 

By using the expansion of the double commutator (2.14), products of commutators and aG 
satisfy the relation, 

[iP^ [P\ aG]]aG = -[V^, aG]^ = aG[V^, [P^, aG]] 

= -AV^VpaGD.DpaG ( 2 . 22 ) 

except both higher order terms and nonessential total derivatives which are not needed in 
our purpose. So the first term in (2.21) cancels with the second term. Hence, hz{x) is 
reduced to 

h2,{x) = 2ig‘^mli:i{x\ ^ — :^VyVp\x)'^^aGDyDpaG . (2.23) 

Here, owing to the algebraic relation (2.5), an extra G^i^ could appear as VyVp = 
\{Vu,Vp} + igG^p. However, since there is the 0{G^) content in the trace, we can ne¬ 
glect the noncommutativity in the matrix element of (2.23). Then the operator Vp can be 
replaced by Pp as we performed in the calculation of H^i^x) before : 
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(2.24) 


{x\ 


{V^ + m. 


_-p p 

2\n' ''' 


A=0 


W = / 


d^P PuPt, 

{^tiY {p^ + mlY 

dfjiu 


2®7r^(n — l)(n — 2)(n — 3)mc*'"' 
Substituting Eqs.(2.19) and (2.24) for (2.23), we obtain 


hYx) = 


ig- 


2 ® ■ 37r^m? 


£ /^a/^b/^c 

Jabc^ ^ ^ • 


Hence, 


H 2 {x) = hi{x) + /i 2 (a;) + /i 3 (a;) 


^9 

2%2 






p /^a/^b/^c 
Jabc^ ^ ^ 


2^ ■ 37r2m2 


Finally, since /i 2 (a^) cancels H^i^x), we get 
77 ( 0 ;) = hi{x) + /i3(a:) 


2%2 


G“G“ 




£ ^a^b/^c 

Jabc^ ^ ^ 


2 ® ■ Sir'^mi 


(2.25) 


(2.26) 


(2.27) 


As expected, the first term in H{x) cancels with the first term in (1.10), which is the 
contribution from noninvariance of the measure and the rest part leads to the divergence of 
the c-quark axial current in the form 


(<9^c^(a;)7^75c(a;)) = - 


tg- 


£ /^a/^b/^c 

Jabc^ ^ ^ 




(2.28) 


We would like to stress that our answer for {( 9 ^cl(a;) 7 ^ 750 ( 0 ;)) is 6 times less than that of 
Halperin and Zhitnitsky [3]. 


III. MATRIX ELEMENTS OF GLUONIC OPERATORS IN EFFECTIVE 

ACTION APPROACH 

We face now the problem of the calculations of varions matrix elements of the gluonic 
operators like those in (1.9) and (2.28). 

First we rederive the Diakonov-Petrov (DP) Effective Action [28] as it was snggested 
in [27]. It is natural to choose the singular gauge for the instantons in describing many 
instanton effects in the propagation of the quarks. In the case of a small packing parameter 
it is possible to do the following sum ansatz for the background instanton field: 

N+ N- 

H ^+) + H A) , (^± = {z±, G±, p±)), (3.1) 

+ 

where Zi, Ui and pi are the position, orientation and size of the Ath instanton. The canonical 
partition function of the iV+instantons and N_ anti-instantons can be schematically written 
as 
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^ N+,N^ 

Zn+,n_ = / detArexp(-\4) d^ZidUidn^pi), 


(3.2) 


where Vg is the instanton-(anti)instanton interaction potential generated by the gluon field 
action and detAr is a quark determinant in the instanton field. The main assumption of 
the instanton model is that Vg is repulsive at small distances between instanton and anti- 
instanton. This should provide the stabilization of the instanton sizes and of the inter- 
instanton distances. We mainly deal with detAr, which describes the influence of light quarks. 

Lee&Bardeen [30] (LB) calculated the quark propagator in a more sophisticated approx¬ 
imation than DP, finding 


detAr = det 5, Bij = imSij + aji, (3.3) 

where atj is the overlapping matrix element of the quark zero-modes generated by 
instantons^. This matrix element is nonzero only between instantons and anti-instantons 
(and vice versa) due to the chiral factor in q ; h e., 

a_+ = -(<h_,oN'9|<h+,o)- (3.4) 

The overlap of the quark zero-modes causes quarks jumping from one instanton to another 
during their propagation. 

Eq. (3.3) implies that for iV+ ^ N_ detAr ~ so the fluctuations of |iV+ — N_\ 

are strongly suppressed due to the presence of light quarks. Therefore we assume = 
N_ = N/2. 

Let us rewrite the detAr following the idea suggested in [31]. First, by introducing the 
Grassmanian (7V+,7V_) vector 


and 


we can rewrite 


D = {ui...un+,Vi...vn_) 
n = {Ui...Un+,Vi...Vn_) 


detAr = 


dfldfl exp (QBQ ), 


(3.5) 


where 


DBG = V,{im + a^)V, = mu+u+ + * E + E('*^+'^-®-+ d" V-U+a+-). (3.6) 

+ - +- 


The product can be rewritten in the form 

'u+rr_a_+ = —(i(9<h_,o'i'-|(*'9)~^|*i9<h+,oh+)- 


(3.7) 


is the solution of the Dirac equation P^±fi = 0 in the instanton(anti-instanton) field 
d.it,/i(3j; C±)- 
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The next step is to introduce iV+, iV_ sources 77 = ( 77 +, 77 .) and iV_, iV+ sources 77 = ( 77 ., 77 +) 
defined as: 


77_ = (i,9<l>_,oV-|,?7+ = {id^+flU+\ 

77 + = | 7 ( 9 <I>+,om+), 77 _ = . 

Then exp(r2a^r2) can be rewritten as 

exp(r2a'^r2) 

= exp J (^—fj{id)~^T]j 

= (^det( 7 ( 9 )) J D'ljjD'ijj'^ exp J dx {x)id'ijj{x) — f]{x)^jJ{x) + ^jJ\x)r|{x)'j (3.8) 

Integrating over Grassmanian variables G and G and taking into account the Nf fiavors 
detAT = 0/ det Bf this provides the fermionized representation of the Lee& Bardeen’s result 
for detAT in the form: 


(3.9) 

X n |n 0"^/+^+[^/>^/]) n 0"^/+| > 

where 

= J {tlj]{x)id^±fi{x]^±)) J d% . (3.10) 

Eq. (3.9) coincides with the ansatz for the fixed N partition function postulated by DP, 
except for the sign in front of V±. Keeping in mind the low density of the instanton media, 
which allows independent averaging over positions and orientations of the instantons, Eq. 
(3.9) leads to the partition function 


det N = J exp yj d^xJ2 

(N+ N- 


Zn = 


exp d‘^xiljHdilj^ w!! , 


where 


W±= f d'^^±Y[(V±[ipj'ipf] + irrif) 

J ^ 




and 


Nr. 


r d*kd^l 


m Nc 


(27r 


(3.11) 


(3.12) 


(3.13) 


The form factor F is related to the zero-mode wave function in momentum space ‘h±(/i:; .C±) 
and is equal to 
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(3.14) 


= [/„(()*•„(«)-/,(«)A-,(«)] , « = iyFp. 


Correlators in the DP effective action 

In quasiclassical (saddle point) approximation any gluon operator receives its main contri¬ 
bution from instanton background. In the following the operator g‘^GG{x) will be considered 
for the illustration of the method. Owing to the low density of the instanton medium, it is 
possible to neglect the overlap of the helds of different instantons. In that case, the matrix 
element of g‘^GG{x) with any other quark operator Q is 


{g‘^GG{x)Q) N = j D'lpD'ip'^ exp d'^xi/jUdi/^j 

X { n + (H'GG+(i)e) + /V_ (W'GG-(i^)e) H-'Cn-'"-'} , 


where 



z) det f 


+i 


m NA 

4^2 p 2 J 


and f 2 {x — z) is dehned as 


(g^GG(x))± = ±f 2 (x - z) 


^ 192p^ 

[p2 + (x - zy]^ 


(3.15) 

(3.16) 

(3.17) 


It is useful to introduce the external held k, 2 { 3 ){x), coupled respectively to g^GG 
{g^fabcG'^&G^). Starting from (3.15) and (3.16), we hnd the partition function Z[k 2 { 3 )\ 
describing mesons[26] in presence of such external held: 


^[<^ 2 ( 3 )] = J exp (-hh[<l>+, $_]), 


(3.18) 


where 

W[^+, <l>_] = J d*x{Wa + Wb- Wc), 


w 


a 


{Nf - 1 ) 


N 

w 


fnM/Met $+1 


+ (<^>+ ^ , 


Wb 


Nc 

4'jl2p2 


Tr{m(<l>+ -I- <h-)}. 


(3.19) 


Wc = X! Tr 
/ 


In 


id + iruf 





2(1 iT's)- 
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The two remarkable formulas 


(ab)^ = J d\ exp(iV/n^ - N + \b) {N » 1) 


and 


exp(A det[iA]) = / (i<hexp 


-{Nf - 1)A (det $) + itr{^A) 


(3.20) 

(3.21) 


have been used here. It is possible to check these formulas by the saddle point approximation 
of the integrals. They were proposed in [29] and we followed this approach. 

The saddle point of the integral (3.18) is located at (*h±)/g = Mf6fg, a self-consistency 
condition for the effective quark mass, i. e., 


4NrV 


d^k 


MjF\e) 


{2ny M]F\k^) + 


= N + 


rufMfVNc 

27r^p2 


(3.22) 


being imposed, which describes also the shift of the effective mass of the quark Mj due to 
current mass ruf. In the following we will neglect ruf, since this model fails to reproduce 
properly the low-energy theorems beyond chiral limit [27]. 

The solution of a self-consistency equation (3.22) in chiral limit correspond to Mq = 
340 MeV assuming the parameters p and i?, which are the values given in (1.11). 

The action IT[<h+,<h_] in (3.19) has imaginary part, in general, which is reduced to 
Wess-Zumino term in long-wave limit {k -C M, where /c is a mesons momentum). 


IV. VIRTUAL CHARM IN B ^ Kp' DECAY 

We apply the Effective Action (3.19) and the formula for the divergence of the c-quark 
axial current(2.28) to the calculation of and compare with the analogous quantity f^\ 

which is dehned in the similar way as in(1.4). These quantities certainly are dehned 
in Minkowski space. The eqs. (1.9) and (2.28) may be easily translated from Euclidean 
to Minkowski space accordingly the above-given prescription. For instance {{G°‘G°')e —^ 
(G“G'“)m) and ((/abcG^G'^G'')^ ^ - (/a6cG“G^G‘')M), which lead to 

i44r> = < oi^(G“G“)„|.,' >. (4,1) 

and 

= - < 0|y3^(/<.i'»G“G'’G')M|i,' > (4.2) 

The phenomenological way of the estimation of the is the application of the QCD-I-QED 
axial anomaly equation together with data on p' 27 decay leads to 

47 ^ = 63.6 MeV, (4.3) 
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which was used in [ 6 ]. On the other hand the calculation of the matrix elements (4.1), (4.2) 
may be reduced to the calculation of the correlators 

< G^G\x) G^G^{y) >, < fa,cG^&G\x) G^G\y) > 

respectively. The calculation of the correlators are naturally performed in Euclidean space. 
These correlators have almost the same dependence on the large relative distances \x — y\ 
and the ratio of these correlators become almost constant, at least in DP chiral quark model 
128 ]. 

In the Effective Action approach, the abovementioned gluonic operators vertices in the 
correlators are generated by differentiation of the Effective Action (3.19) over ^ 2 ( 3 ), which 
leads to the vertices 

* / 2 ( 3 ) NJ^ 75^((4>+ + $-) + ($+ - 4>-)75) (4.4) 

So, the calculations of and also may be reduced to the calculations of the similar two- 
point 75 -singlet correlators with additional form-factor /3(g) or /2(g) (3.17), which are simply 
the momentum representation of the instanton contribution to the operators fabcG°'G^G^ 
and (see Eq. (3.17)), respectively. They are dehned as 

= —1536p® J d^x exp{iqx) {p^ 

/ 2 (g) = 192 d^x exp{iqx) {p^ + x‘^)~^, (4.5) 


and 

/2(0) = 327r2, /3(0) = -^/2(0) 

It is easy to reduce the integrals in (4.5) to the Bessel functions as 


J d'^x expipqx) (p^ -|- x^ 


—n 


{n — 1 )! 



(4.6) 


(4.7) 


With the Effective Action (3.19) we calculate everything in Euclidean space with further 
analytical continuation to the Minkowski region. 

On the other hand the ratio of abovementioned correlators in the Effective Action ap¬ 
proach equal the ratio of form-factors /3(g) and /2(g)- This ratio has a weak dependence on 
the argument (like 10 % on the scale g ~ P~^)- As a result, this provides the possibility to 
calculate the ratio at small Euclidean g^ within this accuracy. 

Now it is clear that with the Effective Action (3.19) the matrix element in (4.2) may be 
reduced to the calculation of the matrix element in (4.1) with an additional factor —So, 
the ratio of Eqs. (4.2) and (4.1) is equal in this model to: 


f(c) 

fW 

Jp 


12 1 
5p2 Qml 


- 0 . 1 , 


(4.8) 
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where the value for p is given in (1.11). By taking into account the estimation (4.3) (we use 

mc{p — nic) — 1.25 GeV on the scale p ~ rric for the numerical estimates), we hnd 

/fo = -6MeV. (4.9) 

This number is close to the one of [5], = 5.8 MeV and the sign and the order of the 

value coincide with the estimations of [6, 8,13]. 

Recently, Shuryak and Zhitnitsky [4] performed direct numerical evaluations of the var¬ 
ious correlators of the operators g^G°'G°', fabcG^G^G^ in the Interacting Instanton Liquid 
Model(IILM). Their calculations lead to: 

= 7GeV^ (4.10) 

(which leads to = 48.3 MeV) and 


mg^fabcG^G^G^W)\ 

mg^G-G^W)\ 


^ (1.5 ~ 2.2) GeV^ 


(4.11) 


The later is somewhat large than their simple estimate for this ratio of matrix elements 
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« (1 ~ 1.5) GeV^ 


(4.12) 


With the use of (4.11) and (2.28) we arrive at 

Ac) 

qL = _o.i7 -0.25. (4.13) 

Au) ^ ^ 

J rj' 

fc) 1 

This ratio gives fy = —8.2 ~ —12.3 MeV at the scale of the size of the instanton p ^ p 
The abovementioned experimental numbers (1.1) are given at the scale p ~ me, which 
is different from the scale of this instanton calculation. The account of the anomalous 
dimension of the g^GGG operator [35] leads to correction [4] 

f^f\p ~ rric) ~ l.bflf{p ~ p~^). (4.14) 

The account of this scale factor leads to 

4f)(/i~me) = -12.3-18.4MeV. (4.15) 

Hence, using (4.11), the result of more sophisticated calculations of Shuryak and Zhitnitsky, 
we get the number (4.15) which is 2-3 times larger than our simple estimation (4.9). 

These numbers (4.9), (4.15) are in agreement with the phenomenological bounds [7,13] 
and almost in agreement in the sign and the value with [6, 8] but six-ten times less than the 
estimations given by [3](see (1.6)). 

By using the numerical analysis of the branching ratio for p'K^ given at [5] (Fig. 17 

of [5] ) we expect that the value of given in (4.15) may provide a more satisfactory 
explanation of the experimental data. 


16 



V. VIRTUAL CHARM IN POLARIZED DIS 


From the definition Eq. (1.7) and from Eqs. (1.9), (2.28) it is easy to find 

NfO^ ~ 

AJ:2mNpi'l5P = {p,s\j^GG\p,s), (5.1) 

• 3 

Ac2mNpi'l5P = - {p, s\fabcG'^G^G''\p, s) (5.2) 

(These formulas are defined in Minkowski space). As it was mentioned before AS = 1 in 
the constituent quark model. On the other hand it was shown that Skyrme soliton model 
gives zero for this quantity. Chiral quark-soliton model (see recent review of this model [38]) 
interpolate between these constituent quark and Skyrme models and give the value of AS 
in the range [36,37] 


AS = 0.3 ~ 0.5. (5.3) 

The calculation of the same quantity by QCD sum rules approach give the same answer (see 
the review [17]). 

Chiral quark-soliton model is essentially based on the Effective Action like (3.19), where 
all the degrees of freedom are frozen except constituent quarks and pions. In that case, by 
restoring of the quark degrees of freedom in (3.19), we find 

ZnM = I Di:^Di:DUexp J (l + 75(«2(3)/))'^'(5.4) 

where 7/'’'® = ^((<l>+-|-<l>_) -|- ($+ —<l>_) 75 ), <l>± = Mexp(±i 0 ) and the usual decomposition 
for the pions 0 = X]i may be used. 

The mass of nucleon is calculated from Euclidean large-distance asymptotes of two-point 
correlator of the composite quark operators T]\[{x) with nucleon quantum numbers and 
certainly by using Zn[0]. 

The nucleon mass receive the contributions from the products of the constituent 
quark propagators in the external pion field and the polarization of the constituent quarks 
vacuum by this field(effective action for the pions) integrated over this field. The saddle- 
point condition in this path integral means topologically nontrivial pion field (ps like famous 
skyrmion (see for example [38]). 

The calculation of the nucleon matrix elements of the any combination of the gluonic 
operators g^G°‘G°' and fabcG°‘G^G^ may be reduced to the differentiation over ^ 2 ( 3 ) of 
the two-point correlator of the composite quark operators T]sf{x) calculated now by using 
^Ar[^2(3)]- 

This is the way to calculate AS and Ac, which are essentially reduced to the nucleon 
matrix elements of the gluonic operators g‘^G°'G°' and g^ fabcG°‘&G'^ respectively. It is clear, 
that these operators lead to the analogous (4.4) vertices 

*/2(3)MF2 7/^®iV-i75. (5.5) 

The form-factors /2 and /s are defined in (4.5). In the present case we need these form-factors 
at = 0 . 
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On the very general ground, it is possible to prove that the calculations of AS with 
using Zj^[k2] leads exactly the result of Wakamatsu and Yoshiki [36]. In that paper authors 
calculated the nucleon matrix element of the operator + d^ 7 ^ 75 (i by using the 

Effective Action like Zn[ 0 ]. 

Let us change variables in Zj^[k2\, Eq. (5.4), accordingly: 

Ip = exp {a{z)'y 5 ) Ip, ip^ = exp (q;( 2 :) 75 ) . (5.6) 

If we choose 

a = 172(^2/2)^^ 

then most important 0{k2) term in Zm[k2\^i\\ appear now in an other form as: 

= J Dtp' J DU 

xexpj dztp^'{z) (^J dx K 2 {x)dz^f,f 2 {x - z)^ ^ tp'{z) (5.7) 

where J is a Jacobian of the transformations of the measure. Such type of Jacobians is 
responsible for the chiral anomaly, in general. In the present case we are dealing with the 
action in (3.19), (5.4), which has a imaginary part, since Dirac operator 

D = id + 

is not hermitian. To calculate the Jacobian J we must properly dehne the measure. 

First of all we must choose the total systems of the eigenfunctions of some hermitian 
operators. It is natural to take the operators: 

D^D = -d‘^ + F^M^ - F‘^M{dW^), DD^ = - + F‘^M{dW^^), (5.8) 

and dehne a set of the eigenfunctions: 

D+D(pn = Xl(pn, DD^^n = (5.9) 

It is easy to show that: 

D (pn ^n^ni 

where \r]n\ = |A„|. 

With this basis we may expand: 

H «n0n, 67®+, (5.10) 

n n 

where and 67 Grassmanian numbers. So, the measure may be dehned as 

Dtp = Y[ dttn, Dtp^ = [Q dbi- 

n n 

The Jacobian J then is given by 

J = exp{J dxa{x)A{x)). (5-11) 
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According to Fujikawa [23], the regularized expression for the anomaly A{x) is given by 


D'^D DD^ 

A{x) = E 0 +(a;)75exp(- —)(f)n{x) + ^ $+(a ;)75 exp(- 


= lim tr 

/i.—^OO 


d'^k 

w? 


exp(—iA;a ;)75 |exp(- 


D+D^ 


+ exp(- 


DD+^ 




exp{ikx). 


(5.12) 


Redehning the variable of the integration in (5.12) k k/fi we arrive at: 

A{x) = tr / exp{-k^)^ 5 [{F^M{dU'^^)f + {F^M{dU'^^+)y] = 0, (5.13) 

J [27rp 


due to the trace over Dirac matrices. So, Jacobian J is equal to one and now it is absolutely 
clear that the calculations with Eqs. (5.7), (5.4) leads to the same result in [36,37]. 

The SU{3) extensions [37] of the model (5.4) will leads to the same result for AS, 
since the valence u, d quarks give the most essential contribution to this quantity, while the 
vacuum quarks are negligible. 

So, the ratio Ac/AS can be easily calculated in the same line as for flf / (4.8) which 

leads to 

Ac/AS = 

With the use of (5.3) it means 


12 


5p2 NfQml 


-0.033. 


(5.14) 


Ac = - (0.01 ~ 0.016). 


(5.15) 


Again the quantity Ac in (5.15) is given at the scale of the size of the instanton. On the 
scale ~ rric the account of the anomalous dimension of the g^GGG operator [35] leads to 
the same correction as in (4.14) and we have 

Ac(p ~ me) ~ 1.5Ac(/i ~ p-^) = - (0.015 ~ 0.024). (5.16) 


VI. CONCLUSION 

The problem of the virtual charm axial current contribution has been reduced to the 
calculations of the specihe gluon operator matrix elements by the application of the operator 
Schwinger method, developed for QCD by the ITEP group. Due to the specihe structure 
of this operator these matrix elements receive the main contribution from the instanton 
background. The DP effective action approach, based on instanton model of QCD vacuum, 
has been rederived and applied to the calculations of such types of matrix elements in chiral 
limit, since the reliable answer may be obtained only in this limit. 

We have calculated the coupling of rj' with virtual charm axial current (4.15). The 
obtained value may provide a satisfactory explanation of the experimental data on the 
branching ratio for p'iP^-decay as was shown recently by Cheng and Tseng [11]. 

They demonstrated an impressive good explanation of these data using our value (4.15) for 

f(c) 

Jri' ■ 
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This approach has been applied also to the problem of the charm content of the nucleon 
spin Ac. We conclude that Ac (5.16) is one order of magnitude smaller than the analogous 
contribution of the strangeness As = —0.11 ± 0.03 [15]. 

We acknowledge the support of the COE program and a partial support by the grant 
INTAS-96-0597, which enables M.M. to stay at RCNP of Osaka University and to perform 
the present study. 
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APPENDIX 


We calculate here the determinant of the Dirac operator which appears from the Gaussian 
integral over c-quarks. This calculation was performed by using of the expansion of the 
effective action in the series of G/m?^ by Vainshtein et ah in Ref. [34], However, we 
show here that careful calculations give a different result from that of Vainshtein et al [34], 
Throughout this appendix, we use the Euclidean convention and consider up to the third 
order of G in the expansion. The contribution of the source which satishes the held 
equation D^G^^ = may be negligible, as being consistent to the main part of this paper. 

As explained in the calculation of the determinant of the Dirac operator, 

det \\V + irricW, must be regularized by the regulator mass M, 


{V + imc){p + iM) 
{p + imc){V + iM) 


(A.l) 


Instead of the direct calculation of V, it is more convenient to calculate the logarithm of 
V. The relation. In det ||A|| = Trln||A||, is satished formally with the inhnite-dimensional 
matrix A, 


\nV 


{V + imc){p + iM) 
{p + imc){V + iM) 


(A.2) 


Note that the symbol Tr in this appendix contains the trace over Lorentz and color indices 
and the integral over the coordinate, i. e. 


Tr(- ■ ■) = tr^+c J d'^xixl ■ ■ ■ |a;) 


(A.3) 


while in the main part of this paper it has been written as Tr = tr^^+c only. 
Let us consider the derivative of InD with respect to the c-quark mass rric, 


1 d 

nic dnic 


\nV 



+ ml + 




+ ml 2 + ml 2 

Iq + h + I3 + ■ ■ ■ , 


+ ml 2 + ml 


+ 


(A.4) 


where the single term vanishes because of the trace of the single 

Calculations of both I 2 and I 3 are performed as we have done in the calculation of 
{ 2 mcC^'y 5 c), except for 75 . Since we need up to 0{G^) in our purpose, each of them is 
calculated as follows: 
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him?) = —Tr 


1 9 ^ i g ^ 9 ^ i 

V^ + ml2 V^ + ml2 V^ + ml2 V^ + m^ 


= -Tr 


{V^ + mlY\2 J 

tr^+c {|^^G(a;)| 


2 ^ ■ ?)Ti‘^mi 


d^xfabcG^G^GYx) , 


(A.5) 


himl) = Tr———^ 

2 + m^2 + mf. 


= T y 


(P 2 + ^2)3 


k)tri+c (o-G)^ + (a;| 


(P 2 + m2)4 


|a;)tr2,+cO'G ■ D^aG 


{V^ + mlY 


VuVYx)tYL+c(^G ■ D^D^aG 


9 

2®7r2m2 


J d‘xG‘G‘(x) + y ^ '^ 2.^1 j d‘xfx^G‘G'-G‘(x) . 


(A.6) 


Here we have used the following relations, 


tiL+c {cjGf = -2H ■ trcGGG = 2^fabcG^G^G^ 
tiL+c {(^Gf = 2hrcGG = 2^G^G^ 

D'^aG = 2iga^^G^aGau 


(A.7) 


under neglecting the source which satisfies D^G^^ = j^. 

The calculation of Jq is rather technical since it needs careful treatments of traces of 
inhnite-dimensional matrices. First, we consider the translational invariance of Jq. It means 
that Iq does not change under a shift of the operator by an arbitrary vector g^. 


himl) = Tr 


p 2 + 77^2 p 2 _|_ ^2 


= Tr 


+ ml + ml 


(A.8) 


= Tr 


{V — qY + ml {p — qY + ml 


Here is the shift operator in momentum space, and is the coordinate operator which 
satishes the algebra (2.1). Now, we expand the right hand side of (A.8) into the series in q^. 
The momentum shift of (P^ + ml)~^ is expanded as 
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^iqX _ 


V"^ + 


^-iqX 


OO 1 

= H rr [iqX , 


n=0 


n\ 


V‘^ + ml 


+ ml 


iqX appears n times 

- 2 ] + 


1 


+ ml 


+ ml 


+ ■■■ 


(A.9) 


When one apply the trace Tr with both sides of this equation, the first term of the right 
hand side is equal to the left hand side and both terms do not depend at all on the vector 
This means that other terms which depend on any power of q are identically equal to 
zero, i. e. 


Tr [ iqX , 


= Tr [ iqX, [ iqX , 


= ... = 0 


+ m 2 J ^ ^ ^ ^ ' V^ +ml 

Consider the coefficient of the term which is the second order of q^: 

1 1 ,, 


for any . (A. 10 ) 


Tr [ iqX, [ iqX , 

= - 2 g^g. 




+ ml + ml 

1 


-4Tr 


{V^ + mlY (jl + mlY 

1 „ 1 


+ m 


:V„ 


lY +ml 


r,, - 


P^lPv 


+ mlY 


(A.ll) 


Since g^ is an arbitrary vector, we may choose the average over the direction of g^, i. e. 

As explained already, the coefficient of g^ is equal to zero. The corresponding 
equation looks as 


Tr 


+ 771^)2 {p^ + mlY 


= Tr 


V, 


{V^ + mlY ^V^ + ml 




p 


(p2 + mlY 


(A.12) 


On the other hand, using the relation (V^ + ml)(V^ + ml) ^ = 1, 

1 1 1 


Tr 


+ ml)‘^ (p 2 + m^ 


= Tr 


1 


2'i2 




P 


+ mlTr 


iV"^ + ml)'^ {p^ + ml)^ 

Hence, substituting it to Eq. (A.12), we obtain the relation 

1 1 1 


{V^ + ml)^ {pi + ml)'^ 


(A.13) 


Tr 


{V^ + ml)'^ (p 2 + ml)"^ 


= ^Tr 


ml {V^ + mlY 




{V^ + m' 




- 1 )^ 


= _Tr _ Z1 _ V 1 "P^ [ -p^ . . . [ -p^ -p 

™2 /''TS2 I ™2^n+3 ' 1 ' ’ L ' ’ L ' ’ ' M. 




(P2 + m2)-+3 ^ ^ 


p 2 

appears n times 
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= ixr 


nil 

1 

^ (p 2 + miy 


+ 


-VJVYVA - 


V^ + ml 


[V\PJ [V\V. 


P^.PA [P\ [P\P,^[ 


(A. 14) 


[PY [P\Pf.]][PY [P^.P^.\ 


(P2 + mlf ' 

To integrate the right hand side over momentum space up to the desired order of G, we need 
only two commutation relations under being neglected, 

[P :Pim] =‘^igPuGuti lA 

[PY[PYP^]] = -^gPaPuD^G,^ - Ag‘^Pe.G^,G,^ + 2tgV^D^G^^, ^ 

since the product of P^ and the higher commutators which could give the 0{GY terms can 
be replaced by the product of these commutators (A. 15) with a suitable sign in the trace 
Tr. Thus, 

1 1 1 


Tr 


+ miy + ml)^ 


g 


Id^xG^G^ix) - / d*xUcG^G’>G^{x) . (A.16) 


2 ^ ■ J 2 ^ ■ 3 ■ 57r^m® 

The integration of this result over squared masses twice gives 

1 1 


roo poo 

= 2 / dm\ / dm^Tr 

.! YYl"^ 


Jm- 

.2 


iP"^ + mlY {p^ + mlY 


g 


2 ^ ■ STr'^ml 
Finally, we obtain the determinant 

rM^ 


J d^xG^G^ix) - ^ 3 . 32 ^ 3 ^ 2^4 / d^xUcG^G^G%x) . (A.17) 


1 ( \ 
T> = exp- / dmP + h^ml) + h{ml)\ 

= exp [ 


g 


2 ^ ■ 3n^ ml 


In ^G“G“(a;) - 


23g^ 


2 '^ ■ 3^ ■ 


U,G^G^G%x) 


(A.18) 


To compare with the result of Vainshtein et ah [34], we needs to obtain the effective 
action written in the Minkowski metric. In our case, the Euclidean effective action is dehned 
as 


SYsf = \riV. (A. 19) 

According to the previous convention for the transition between the Minkowskian and the 
Euclidean space, 

{GG)e = {GG)m , {GGG)e = -{GGG)m • (A. 20 ) 

Then the Minkowskian effective action is 

Sih =Jd‘x {-^ lu ^trcGG(x) + 3 f£.„. trcGGG(x)| . (A.21) 

Here, tr^ is explicitly written in order to compare with the convention of Vainshtein et ah 
[34], The coefficient of the second term is different from their result [34] with factor 23/2. 
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